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ABSTRACT 


Under the arid conditions prevailing at the end of the dry season (near 
the Autumn equinox) in the western Negev/ northern Sinai region, vegetation 
causes a sharp increase relative to bare soil in the daytime sensible heat 
flux from the surface to the atmosphere. Two distinct mechanisms are 
involved: the increase in the surface abosorptivity (co-albedo) and a 
decrease (as compared to bare soil with the same absorption of insolation) in 
the surface (soil/plants) heat flux. Based on the values of the soil albedo 
and of the plant density in a fenced area (an exclosure) in the Sinai, the 
increase in the surface absorptivity by vegetation fully protected from 
grazing is by a factor E da of 1.23 at solar zenith angle of 30° and by a 
factor of 1.33 at solar zenith angle of 60°. Neglecting the latent heat 
flux, which corresponds to less than 0.5 mm of water evaporated per day, this 
mechanism increases the sensible heat flux in the same ratios. These desert- 
fringe plants, mainly thin, predominantly vertical grey/brown twigs and 
stalks, readily transfer the absorbed solar radiation to the atmosphere. The 
intimate aerodynamic contact between the plants and the air flow can be 
inferred from the fact that the plant temperatures are relatively close to the 
air temperatures. The heat storage in the plants is considered negligible. 
The soil heat flux is significantly reduced when compared to that into bare 
soil with the same surface albedo as the plants/soil surface. Assuming that 
for the bare soil the daytime soil heat flux is half of the net radiation 
(from measurements by Clothier et al., 1986), this soil heat flux reduction 
translates according to our model into a daytime enhancement E^ of the 
sensible heat flux by a factor of 1.27 at solar zenith angle of 30°, and by a 
factor of 1.47 at solar zenith angle of 60°. The combined effect of the 
plants, E^, by both the albedo and the soil heat reduction can thus 
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increase the sensible heat flux by a factor of 1.56 at solar zenith angle of 
30° and by a factor of 1.95 at solar zenith angle of 60°. These calculations 
are for the direct solar beam only; a less steep increase with solar zenith 
angle occurs when the scattered radiation is considered. 

By increasing the sensible heat flux to the atmosphere through the albedo 
and the soil heat flux reductions, the desert-fringe vegetation increases the 
daytime convection and the growth of planetary boundary layer. When 
moisture moves from outside the dry region, probabilities of precipitation are 
then enhanced. The effects can be expected to be especially large during the 
morning and afternoon hours, that is, at moderately large solar zenith angles. 
Conversely, removal of vegetation by overgrazing, by reducing the sensible 
heat flux, tends to reduce daytime convective precipitation, producing higher 
probabilities of drought conditions. 

This assessment of overgrazing (or conversely, of protecting the 
vegetation over a large region) is based on observations in the Sinai/Negev, 
where the soil albedo is high (about 0.42) and where overgrazing produces an 
essentially bare soil. The calculations should not be automatically applied 
to all arid regions. The effects would be smaller when calculated for soils 
with a lower albedo. In the Kalahari and in northern Mexico, bushes or trees 
remain in the overgrazed areas, and strong climatic effects of overgrazing 
cannot be postulated. Even if our assessment for Sinai/Negev does not 
quantitatively apply throughout Africa, the current practice in many African 
countries of maintaining a large population of grazing animals, can contribute 
through the mesoscale mechanisms described here to reduced daytime convective 
precipitation, perpetuating higher probabilities of drought. Time-of-day 
analysis of precipitation in Africa appears worthwhile, to better assess the 
role of the surface conditions in contributing to drought. 
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1 . Introduction 


During the last fifteen years, vast regions of arid and sub-arid Africa 
have experienced incidents of severe drought. While drought in these regions 
had been recorded in the past (Nicholson, 1979), the recent persistence of 
drought is puzzling. The statistics of precipitation are such that the failure 
of only one or two intense rainstorms to materialize in a given year signifies 
a drought year (S. Nicholson, personal communication). 

Addressing this problem area, we discuss here how denuding the surface of 
its natural vegetation can reduce the convection and thus reduce the 
probabilities of convective precipitation. This study is aimed to be 
applicable to the conditions prevailing at the end of the dry season 
(September/October) in northern Sinai/toestern Negev. The analysis is based 
on measurements by Idso et al. (1975) and Clothier et al. (1986) of the net 
radiation and the soil heat flux; on a simple model of absorption of solar 
radiation by the plants/soil surface, developed specifically for the desert- 
scrub growing on the extremely bright soil of the northern Sinai (Otterman 
1981a); and on interpretation of ground measurements and satellite data over 
an exclosure in this region (Otterman 1981b; Otterman and Robinove, 1982; 
Otterman and Tucker, 1985). 

The difference between the net radiation balance R_ and the soil heat 
flux G, measured and discussed by Idso et al. (1975), constitutes the primary 
term for evaluating the potential evaporation from which the actual 
evaporation (the latent heat flux from the surface to the atmosphere) can be 
assessed under a set of assumptions or measured conditions. In our analysis 
we apply the data obtained by Idso et al. (1975) to a different purpose, to 
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assess the sensible heat flux from a dry surface. It is thus appropriate to 
discuss first the climatological circumstances that prevail in Sinai/Negev at 
the Autumn equinox, before the first rain, which justify this approach. 

The Negev and the northern Sinai lie at the fringe of the great desert 
belt. Based on changes in the distribution pattern of various woody species 
during the last 5000 years, the climate of this region has apparently become 
drier and warmer (Waisel, 1986). The rains in Sinai/Negev come predominantly 
during the four winter months, November - February, with precipitation ranging 
from 100 mm y~^ in the northern Sinai to 200 mm y~* in the western Negev (to 
the north of the northern Sinai). Dew in the summer possibly can contribute 
20 mm y - ^ of additional moisture. Considering strong run-off, no more than 
120 mm is available for evaporation/evapotranspi ration during the eight dry 
months March - October, that is, no more that 0.5 mm/day as the average. 
Toward the end of the dry season the evaporation should be less than this 
average figure. Analyzing in our study the differences between vegetated and 
denuded areas, we focus on the ratio of the sensible heat flux from a 
vegetated surface to that from a bare soil. The latent heat flux from either 
surface is neglected. The calculated ratios apply more accurately to the 
combined sensible and latent heat fluxes. If we assumed the evaporation rate 
for each surface as 0.5 mm day - ^, the enhancement ratios for the sensible heat 
flux by the plants would be slightly higher than the ratios reported here. 

Our aim in this study is to treat in a realistic manner the absorption of 
the solar radiation by a vegetated surface in a desert-fringe region and the 
transfer to the atmosphere of the sensible heat. Daylong measurements of the 
radiation temperature of desert scrub and of soil in the Sinai exclosure 
(Otterman, 1981a) clearly indicate that the temperatures of these two surface 
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components as a function of time of day follow distinct, though similar, 
patterns. Both components show radiation temperatures rising towards noon and 
declining in the afternoon. However, these plants, predominantly dry twigs 
and stems protruding approximately vertically from the bright soil (for a 
photograph of the exclosure, see Otterman and Tucker, 1985) intercept most of 
the direct solar beam when the sun is near the horizon, and therefore in the 
early morning warm up more rapidly than the soil. Later in the morning, the 
soil warms up more steeply, at which time the plant radiation temperatures 
become lower than those of the soil, closer to the air temperatures (Otterman, 
1981a). 

These measurements suggest that we analyze the absorption of solar 
radiation at the surface separately by the plants and by the soil - and treat 
separately the transfer of the absorbed radiation from these two surface 
components to the atmosphere or to the soil. Such calculations are presented 
here, based on a simple model previously developed by Otterman (1981a). 
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2. Absorption of solar radiation by the plants and by the soil 

In this section we analyze how the combined absorption at the surface is 
increased by the protruding plants as a function of the time of day. In 
reading this and the following section, one can visualize the plants as 
floating at a height of several meters, rather than attached to the soil, and 
thus -one can regard the absorption by the plants as equivalent to that by a 
chaff layer. The model (Otterman, 1981a) represents the plants as thin 
vertical cylinders and the soil as a horizontal plane reflecting according to 
the Lambert law. The reflectance of the plant components is by a factor of 
more them 3 lower than that of the Sinai soil in the vicinity of the exclosure 
(Otterman, 1981b). We therefore make here a further simplifying assumption 
(as compared to the model presented in Otterman, 1981a) that the plants are 
total absorbers, of zero reflectance and zero transmittance. Based on these 
assumptions, the absorptivity b dg of the soil for the direct solar beam is: 

b dg (s ' r i' 9 o ) = (1_r i ) ex P( _s tane Q ), (1) 

where © Q is the solar zenith angle, r^ is the Lambertian reflectance of the 
soil (i stands for interstices) and s is the projection on a vertical plane of 
plants per unit area. Interception of solar radiation by the tree trunks was 
formulated in a similar way by Federer (1968) in his two-level forest canopy 
model. Our protrusion parameter s is identical to the silhouette parameter 
introduced by Lettau (1969) into formulation of surface roughness (Lettau 
denoted his silhouette parameter as s/S , where s is the silhouette of an 
average obstacle and S is the horizontal area it occupies). In Eq.(l), the 
term exp(-s tan© Q ) denotes the fraction of the direct beam that penetrates 
through the plants and reaches the soil and thus 1- exp(-s tan© Q ) is the 
fraction intercepted and absorbed by the plants. The plants absorptivity % 
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for the direct solar beam is: 


b dp * s ' r i' ® 0 ) =1-ex P( -s tan© 0 )+I(s) exp(-s tan© Q ), 


( 2 ) 


where the last term describes the absorption by the plants of the fractional 
flux r^ exp(-s tan© Q ) reflected from the soil. I(s) is the fraction of 
that flux intercepted and absorbed by the plants: 


fft / 2 - 

I(s) = 2J cos© sin© ^l-exp(-s tan0)Jd©, 


(3) 


where © is the zenith angle of a ray reflected from the soil. 


The absorptivity of the surface b d for the direct beam (the co-albedo) is 
given by the sum of the plants absorptivity and the soil absorptivity: 

b d (s, r^ © o )=l-r i |l-I(s)^ exp(-s tan© 0 ) (4) 

The absorptivities b dg , b dp and b d , and the ratio P d of the plant 
absorptivity to the combined surface absorptivity 

P d (s ' r i' 9 o )=b d/ (b dp +b dg ) (5) 

are now calculated, using 0.42 for the soil albedo r^ This value for r i is 
based on the measurements in the northern Sinai and at the Sinai/Negev 
boundary from Landsat (Otterman and Robinove, 1982), and from the ground 
(Otterman, 1981b). For the protrusion parameter s we adopt the value of 0.2, 
as inferred for the 1974 Sinai exclosure from Landsat observations by Otterman 
and Robinove (1982) one and two years after the fencing off. The same value 
of s was inferred by Otterman and Tucker (1985) from NOAA AVHRR data seven 
years after the fencing off. With these values of r i and s, we tabulate 
b dg' b d p ' b d and p d in the top four rows of Table 1 (below © q ) for solar 
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zenith angles 30° to 70°. 

We note that the plant absorptivity increases steeply with increasing 
© Q , while b^g decreases. Their sum, b^, increases with 9- o . The fraction 
increases steeply with 0 Q . 

In the fifth row of Table 1 we tabulate the ratio E da of the plants/soil 
surface absorptivity b^ to that of the same soil bare of vegetation, that is, 
to the co-albedo of l-r^=0.58. The ratios b^/0.58 increase from 1.23 at 30° 
to 1.33 at 60° and 1.40 at 70°. In these calculations solely for the direct 
solar beam, the dependence on the solar zenith angle O q is somewhat 

exaggerated, compared to the situation in which the scattered radiation is 

considered. The advantage of Eq. (1), (2) and (4) is their simplicity as 

compared to the equations for the global radiation (direct and diffuse) 

presented later. The vertical plants would enhance the surface absorptivity 
by somewhat higher ratios, if the anisotropy of the reflection from the soil 
( reflection at low elevation angles larger than that according to the Lambert 
law) is considered (Otterman et al., 1987). 

A simplistic representation is made in this section that the nature of 
the plants/soil surface is (except for the lower albedo) identical to the bare 
soil surface. This would be the case if the dark plants would be embedded in 
the top layer of the soil, thus reducing the albedo without forming a non- 
planar structure. Under this representation, the difference R n -G between the 
net radiation and the soil heat flux can have the same proportionality factor 
to S n , the solar radiation absorbed at the surface, in the case of the 
vegetated (low-albedo) surface as in the case of the bare (high-albedo) 
surface. Idso et al. (1975), for instance, report the difference R^-G as a 
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fraction of S n * On the premise that for given wind conditions and thermal 
emission from the atmosphere R n ~G is a specified fraction of S n> the bj/(l-r^) 
absorptivity (co-albedo) ratio represents the enhancement of the sensible (and 
latent, to the extent moisture is available) heat flux to the atmosphere. 

We now repeat our calculations taking into account the scattered 

(diffuse) solar radiation. Because our surface model is azimuth independent, 

we have to consider only that the scattered radiation illuminates the surface 

generally at a different effective zenith angle them 0 Q . This effective 

zenith angle, denoted 9 , can be expected to be about 50°. It can be assumed 

that © x does not change when 0 Q changes, if scattering by the aerosols into 

aureole region are excluded (The intense forward scattering by aerosols into 

the narrow cone around the direct beam can be simply added to the direct 

beam) . The absorptivity b by the soil of the scattered radiation is thus 

xg 

b xg ^ s ' r i' V = ( 1-1 V ex P( -s tan0 x ) (6) 

and that by the plants 

b xp (s ' r i' e x )asl- f“ I(S,r iJ exp(-s tan© x ) (7) 

To assess the relative magnitude of the solar radiation absorbed at 
different solar zenith angles and different times of the day, we analyze the 
absorption rates of the solar radiation by the soil, B^, and by the plants. Bp 
(rather than the respective absorptivities) . The solar radiation outside the 
atmosphere is given as cos0 Q . At the surface, the direct solar beam f^ is 
given as 

f d (0 o' r s' T a ) = cos0 o expf-(T s+ T a )/cos0 o ] , (8) 
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where T g is the scattering optical thickness and is the absorption 
optical thickness. The scattered irradiance flux f is: 


( 9 ) 


f x (e o' T s ,T a )=COS0 o ex P(- T a / cos0 o ) fl-exp(-T s /cos© o )J/2, 


where the simplified single-scattering approach is used, which is a 
satisfactory approximation under low optical thickness conditions (Otterman, 
1979), and where the intense scattering by aerosols into a narrow cone around 
the direct beam is added to the direct irradiation (scattering outside this 
cone is assumed to be equally divided between directions above and below the 
horizontal plane). The soil absorption rate is: 


V^i-V - f d b dg + E X b xg 


(10) 


and similarly the plant absorption rate B is: 

r 


- f d b dp + £ x b xp 


( 11 ) 


Using the values T a =0.1 and T s =0 . 1 , we compute B^, B^, the ratio P of 
the plants absorption rate to the surface absorption rate: 

B_ 


P(s,r i ,© o ) 


B + B 
P 9 


(12) 


and the enhancement E of absorption by the plants/soil surface to that by the 
bare soil: 


B + B 

£ 3_ 


(f d+ f x ) U-r.) 


(13) 


These four functions are tabulated in the top rows of Table 2 (below 0 Q ). We 

note first that at 0 of 70°, B is more than two thirds of B at 30°, whereas 

° P P 

B at 70° is only one fifth of B at 30°. The surface absorption, the sum B„ 

<3 9 g 
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+ Bp, is at 9 q - 70° about one third of the surface absorption at 30°. 
Comparing now P and in Table 2 to and E^ a in Table 1, we note that these 
ratios of absorption rates (or of absorptivities) are the same at © Q of 50°, 
which is to be expected, inasmuch as we assumed © x of 50°. At solar zenith 
angles smaller than 50°, the ratios P and are higher than the ratios P^ and 
E da for the direct beam only (the differences between the ratios E da and E q 
for 0 Q of 30, 40 and 45° show up only in the fourth decimal place); and 
conversely are lower when Q q is larger than 50°. The enhancement E a 

at 70° stands now as 1.38, instead of 1.40 for the direct beam. In view of 
the many approximations used in these calculations, these numerical results 
suggest that under low levels of aerosols the calculations for the direct beam 
only are a satisfactory simpler substitute for the more accurate but more 
involved calculations for the global radiation. 

The same functions are plotted in Fig. 1 vs. time of day t in hours (t=6 
at sunrise, t=12 at noon). At equinox, the solar zenith angle © Q is given as: 

cos© o = cos|flU2-t)/12|cos L, (14) 

where L is the latitude. L is 31.5° at the Sinai/Negev boundary near the 
Mediterranean coast. The solar zenith angle, plotted in Fig. 1 (at the top), 
decreases from sunrise to the value of 31.5° at noon. Both absorption rates 
increase in the morning but at different slopes: Bp is larger than B^ until 

about 0800h, but at noon B^ is larger than Bp by about 2.5 (see the middle 
graphs). The ratios E. and P decrease monotonically to noon (see the lower 

a 

graphs ) . 
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3. The second role of plants: the enhancement of the sensible heat flux to the 
atmosphere through the soil heat flux reduction. 

Clothier et al. (1986) report that vegetation sharply reduces the soil 
heat flux. This reduction constitutes the second role of the plants: the 

heat not consumed by the day-to-night storage in the soil becomes available as 
a daytime heat input to the atmosphere. Assuming that the storage of heat in 
the dry, thin plants is negligible, we apply our model to formulate the 
reduced soil heat flux and the enhanced transfer of the absorbed solar 
radiation from the surface to the atmosphere as a function of time of day. 

In the previous section we assumed that the R n ~G difference between the 

net radiation R n and the soil heat flux G has the same proportionality factor 

to the solar radiation S n absorbed at the surface for both types of surface 

analyzed here. Thus, the ratio E_ of the differences R -G for a low-albedo 

ct n 

surface to a high-albedo surface is equal to the ratio of the respective 
absorptivities (co-albedo). We regard this assumption as simplistic because 
our vegetated terrain should be analyzed as a complex (nonplanar) surface. It 
can be reasonably assumed that the thin nearly vertical twigs and stalks are 
so thoroughly exposed to the air flow that the solar radiation they absorb is 
immediately released (as sensible heat flux) to the atmosphere. The intimate 
aerodynamic contact between the plants and the air flow is indicated by 
relatively small differences (when compared to the soil - air temperature 
differences) between the plants' radiation temperatures and the air 

temperatures, as reported in daylong measurements in the Sinai exclosure 
(Otterman, 1981a). 

Looking at these effects from a different viewpoint, the plants shadow 
the soil and thus reduce the solar radiation absorbed by the soil. The 
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approach in this section to evaluating R fi and G is based on the assumption 
that . the net radiation R n has the same proportionality factor to the solar 
radiation absorbed by the surface (that is, by the soil and by the plants, if 
any) in both the low-albedo (plants/soil ) and the high-albedo (bare soil) 
case, whereas a different relation holds for G. The longwave radiation losses 
from the surface to the atmosphere are thus assumed to have the same 
proportionality factor to the absorbed solar radiation in the plants/soil case 
as in the bare soil case. In making this assumption, we err on the 
"conservative" side (we underestimate the sensible heat flux from the 
plants/soil surface), because actually the soil temperatures are reduced in 
the plants/soil surface as compared to the bare soil case (with the same solar 
radiation S n absorbed at the surface ) and the plant temperatures are most of 
the day lower than these (reduced) soil temperatures. Moreover, the air 
temperatures over vegetated areas would tend to be higher due to reduced soil 
heat flux (see below), which would tend to reduce further the longwave 
radiation losses. 

The soil heat flux G is assumed to have the same proportionality factor 
in both cases to the solar radiation S r absorbed solely by the soil. S r 
equals S n in the case of bare soil, but in the case of plants/soil, S r is only 
a fraction of S n , given by the ratio of the soil absorption rate to the total 
surface absorption rate: 



The symbols without the subscript d indicate that the scattered solar 
radiation is considered. In terms of the absorptivities for the direct beam 
only, the expressions are simple: 


14 



( 16 ) 


(1-r^) exp(-s tan© 0 ) 

1-r. [l-I(s)] 

In this formulation the influence of thermal emission from above on the soil 
temperatures, and thus on the soil heat flux, is not considered. This 
emission originates in part from the plants rather than from the atmosphere, 
(predominantly so at large viewing zenith angles). The plant temperatures are 
approximately equal to the air temperatures near the surface. At any usual air 
temperature profile this high emission would tend to increase G as compared to 
the case of bare soil. The effect is small however, inasmuch as at large 
zenith angles the atmospheric emission originates anyway effectively from 
near-surface levels. 


S dr _ b dg 
S dn b dg +b dp 


Clothier et al. (1986) relate through a series of measurements under 
various soil moisture conditions the relationship of G to R n . Their findings 
can be represented for our purpose as an approximate relationship: 


G 


9 0 R n 


(17) 


with the values of g Q for dry soils of 0.4 or 0.5. Similar measurements were, 
conducted in the Negev by Fuchs and Hadas ( 1972 ) . 


Let Eg denote the enhancement by the plants of the sensible heat flux 
solely thorough the soil heat reduction (in addition to the enhancement E_ 

cl 

through the albedo reduction). We can state that: 


V^i'VV 




R -G 
n 
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p 


(18) 


HV 11 , : B p> B q + V (1 -9p> 1 + Jo_ 

1 - % B g + B p 

where g r is the ratio of the reduced soil heat flux G r in the plants/soil 

surface to G, the heat flux into the bare soil for the same absorption S n of 

solar radiation at the surface . These statements for indicate that the 

absorption of solar radiation by the plants B^ is by a factor 1/(1 - g Q ) 

(that is, by a factor of 1.67 and 2.0 for g Q of 0.4 and 0.5 respectively) more 

effective in producing sensible heat flux than the absorption by the soil. 

In Fig. 1, a scale is provided on the right, to assess E from P when g =0.5. 

y u 

The enhancement ratio E^ for the direct beam only (subscript d) is 

tabulated for g =0.5 and g = 0.4 in rows 5 and 6 of Table 1. The ratio E 
3 o 3 o g 

for the global (direct and scattered) radiation is tabulated in the same rows 

of Table 2. We note that the enhancement ratio E for g =0.5 is 

substantially higher than the enhancement E & by the reduction of the surface 

albedo. The combined enhancement ratios, the products E^E^g (for the direct 

beam) and E E (for the global radiation), are tabulated in the last two rows 
a g 

of Table 1 and Table 2, respectively. 

Inasmuch as the second role of the plants, of reducing the day to night 
storage of heat at the surface, produces daytime effects higher than those 
solely by the albedo reduction (if g Q is 0.5), the combined mesoscale effects 
can be more than double those hypothesized before (Otterman, 1974) on the 
basis of the albedo reduction alone. 

There is an additional impli cation, for the need to support the albedo 
measurements by the surface temperature measurements. The transport of 
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sensible heat from the surface to the atmosphere is proportional to the 
product of the temperature gradient above the surface by the turbulent 
transfer coefficient for heat, see for instance Eq. (11.9) in Hillel (1980). 
If it is assumed that the aerodynamic contact of the low-albedo surface with 
the air flow is the same (the same turbulent transfer coefficient for heat) as 
that of the high-albedo surface, the increased sensible heat transfer to the 
atmosphere for the low-albedo surface should be evidenced by higher surface- 
to-air temperature differences (if both surfaces are dry). That means higher 
surface temperatures if air temperatures are the same. When early Landsat 
images showed in the spectral bands of the Multispectral Scanner System a high 
spectral albedo of the overgrazed Sinai vs a much darker Negev, this 
consideration prompted an instrumented aircraft flight back and forth across 
the Sinai/Negev boundary. The August 31, 1973, 1400 h airborne measurements 
in the broad 8-11 pm thermal band showed the radiation temperatures of the 
vegetated Negev to be 4 to 5 K higher than those of the overgrazed Sinai (in 
the area of the highest contrast, where there was a fence at the boundary). 
Because of the narrow field of view of the PRT-5 radiometer and the vertical 
structure of the vegetation, in the Negev the instrument views predominantly 
the interstices of soil between the plants. (Our model, which predicts that 
plants are not visible when viewing exactly from zenith, should not be applied 
here). From these measurements, the rainfall reduction mesoscale effects of 
overgrazing were formulated (Otterman, 1974) to parallel the study of Black 
and Tarmy (1963), who suggested tarring the Mediterranean coastal areas to 
increase the surface daytime temperatures and thus increase convection and 
precipitation. 

The interpretation of these radiation temperature measurements in the 
broad thermal band, which included the restrahlen band, was faulted for not 
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taking into account the reduced emissivity of sand in the restrahlen band 
(Jackson and Idso, 1975). Measurements of radiation temperatures at 11.5 pm, 
outside the restrahlen band^from NQAA AVHRR (early morning passes) always 
showed higher radiation temperatures in the vegetated exclosure than in the 
surrounding overgrazed terrain. The differences ranged from 0.1 K to 2.9 K. 
No correlation could be found with wind speed or time from sunrise, to explain 
the spread in these temperature differences (Otterman and Tucker, 1985). 
Thus, the temperature measurements do suggest a higher surface-to-air 
gradients over desert-fringe vegetated area than over adjacent bare soil. 
However, these higher temperatures and higher gradients are not required as 
evidence of the higher sensible heat flow from the surface if the plants 
facilitate the heat transfer to the air flow as compared to that from the bare 
soil. This facilitated transfer can be readily inferred by considering that 
in our model the area of contact with the airflow of plants per unit surface 
area is # s. Thus, in the plants/soil case with s=0.2, the combined area of 
contact (of soil and of plants) is larger by a factor 1.63 than in the case of 
the planar surface of soil without plants. 
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4. Discussion and conclusions 


Our analysis suggests that reduction of the soil heat flux can sharply 
enhance the daytime sensible heat flux to the atmosphere. The enhancement is 
especially strong at large solar zenith angles. For the high soil albedo 
measured in the Sinai/Negev from the ground (Otterman, 1981b) and from Landsat 
(Otterman and Robinove, 1982), and the high value of the soil heat flux, G = 
0.5 R n , taken from the measurements of Clothier et al. (1986), the enhancement 
is by a factor of 1.6 at large solar zenith angles and by 1.3 at moderate 
zenith angles. It is thus significantly larger than the enhancement by the 
surface albedo reduction. The time-of-day representation of the reduced soil 

heat flux adopted here has not been demonstrated experimentally it is only a 

plausible application of our model (Otterman, 1981a). The substantial 
magnitude of the reduction of the soil heat flux by plants has been firmly 
established by measurements of Clothier et al. (1987). They report G/R n 
ratios in the 0.4 to 0.5 range for bare dry soil, and about 0.15 for dry 
vegetated terrain. These results can be translated (without using our model) 
into enhancement of the R n ~G differences by factors of 1.4 to 1.7. Even the 
highest ratios that we present (at solar zenith angle of 70°, see Table 2), 
are lower than that. It would be very interesting to conduct the measurements 
of soil heat flux for vegetated and bare surfaces in desert-fringe regions. 

The combined effects of the albedo and soil heat flux reduction are about 
2 at large zenith angles. Such effects can be expected to cause a significant 
increase in the daytime growth of the planetary boundary layer ( PBL ) . A study 
by Berkofsky (1977), in which the vertical velocities at the top of the 
boundary layer are related to the albedo and the sensible heat flux, is 
pertinent to the problem area. We do not directly address the PBL growth, 
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which should be studied by daylong simulation of the surface-atmosphere 
interactions. Such a simulation is currently being attempted (Otterman and 
McCumber, in preparation) . 

We assume the R n /S n ratio to be the same for the plants/soil surface as 
for the bare soil. In reality the ratio for the plants/soil should be 
appreciably higher, because the longwave loss from the surface to the 
atmosphere is lower: in the plants/soil surface the soil temperatures are 
reduced as compared to bare soil case (for the same S n ) and the plant 
temperatures are most of the day significantly lower (closer to the air 
temperatures) than these reduced soil temperatures. This assumption leads to 
underestimating the enhancement of the sensible heat flux by the plants. 
However, it should be recognized that a significant fraction of the longwave 
loss from the surface is absorbed within the boundary layer, and this heating 
can produce the same effects as the sensible heat flux. 

We do not analyze here how the atmospheric temperatures rise as a 
consequence of the higher sensible heat flux due to the reduced soil heat 
flux. A rise in the atmospheric temperatures would change the thermal 
radiation exchange between the surface and the atmosphere, and this in turn 
would affect the soil heat flux and the sensible heat flux. Without 
considering such interactions, the enhancement ratios that we calculate should 
be considered as only preliminary, with fuller evaluation to be obtained in 
daylong simulation. 

The difficult question is, how the enhancement of the daytime sensible 
heat flux affects the cloud formation and the probabilities of convective 
precipitation. Our discussion is certainly not quantitative in terms of 
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rainfall. The analysis of Black and Tarmy (1963), in which elevated surface 
temperatures of a tarred (low-albedo) surface are related to a large increase 
in the precipitation, might not be quantitatively appropriate (Otterman, 
1974). A careful statement by Malkus (Joanne Simpson) is appropriate here: 
"there exists at least indirect evidence that useful water might be obtained 
in some localities from increased cumulus showers in weakly disturbed or 
undisturbed conditions where skies are generally sunny enough for a thermal 
mountain to be effective" (Malkus, 1963). 

Our analysis deals with mesoscale phenomena, and we do not address the 
large scale effects of regional high albedo, that produce descending 
airmasses over the desert, as analyzed by Charney (1975). Such descending 
ai masses and the associated inversions are important characteristics of the 
climate in Israel and the Sinai throughout the summer, see for instance Shaia 
and Jaffe (1976). The albedo effects accrue in the daytime only, but can be 
integrated over a day or over several days. The soil heat reduction effects 
accrue also in the daytime. However, for this effect the integration over 24 
hours produces practically a null difference in the sensible heat flux. In 
other words, if heating from the plants/soil surface does not produce squall 
lines by evening, the effect of the reduced soil heat flux is null because the 
reduced heat flux from such a surface at nighttime (as compared to that from 
bare soil) would tend to produce at sunrise an atmospheric profile less 
conductive to cloud formation. (Likely non-linearities of the boundary layer 
growth in terms of the sensible heat flux make this statement only 
approximately true ) . 

The effects of soil heat flux reduction on convection and convective 
precipitation is limited to daytime, and can be expected to be especially 
strong in the morning and afternoon hours, at large solar zenith angles. The 
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precipitation in the arid regions is indeed a predominantly intense convective 
daytime phenomenon and this predominance increases with increased aridity 
(Sharon, 1972; Kutiel, 1977). In the western and northern Negev the day/night 
partitioning is somewhat complicated: the early rains, October - November, 
show strong daytime predominance, while the opposite is true for the later 
rains (Otterman and Sharon, 1979). The interesting question in the case of 
the African drought is, at what time of day should the rains that failed to 
materialize in a drought year, have occurred? A study of the diurnal 
variations in convective activity and precipitation during GATE (McGarry and 
Reed, 1978) indicates that in the sub-Sahara region much of the rainfall is 
produced by squall lines that typically form in the afternoon as a result of 
the daily heating cycle. 

Our assessment of the effects of vegetation (or the consequences of its 
removal) is believed to apply quantitatively to the Sinai/Negev region, where 
the soil albedo is high, where overgrazing results in an essentially bare 
soil, and where outside of the desert bloom period the vegetation contains 
very little moisture. These dry twigs and stalks are more like stubble them 
like the green vegetation in more humid regions. Our calculations should not 
be quantitatively applied to all arid regions. In the Kalahari and in northern 
Mexico, bushes or trees remain in the overgrazed area. The effects of 
overgrazing must be much weaker there, both in terms of the albedo and of the 
soil heat reduction. The trees and bushes in such less arid (compared to 
Sinai/Negev) areas do contain appreciable moisture, and the daytime heat 
storage in these plants cannot be regarded as negligible. Finally, the 
effects would be smaller in a region where soils have an albedo lower than the 
0.42 found in Sinai/Negev. 
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The Sinai/Negev vegetation is not unique, however. The desert shrub 
rangelands in Iraq exhibit the same vertical structure in clumps of thin twigs 
and stalks, dark against the background of bright soil, as shown in numerous 
photographs presented by Thalen (1979). Thalen also reports the rapid and 
dramatic recovery of vegetation in overgrazed land, simply by forming an 
exclosure. This confirms our findings (Otterman, 1981b; Otterman and 
Robinove, 1982) and our opinion, that the vegetative cover (outside the desert 
bloom period, during the rainy season) is primarily a function of the grazing 
practices, and not of past precipitation. 

It is thus an open question to what extent our analysis of possible 
mesoscale effects of overgrazing applies to the problem of African drought. 
Applicability of our calculations can be advocated most strongly for the 
extremely arid regions. Under conditions of extreme aridity, the "mean" 
rainfall does not convey the same meaning as in more humid regions, because of 
the large year-to-year fluctuations. We should regard precipitation in terms 
of probabilities of zero, one, two or three rain events in a year. In such 
regions, the current practice in many African countries of maintaining a large 
population of grazing animals results in reduced vegetation density. The 
reduced vegetation, through the mesoscale mechanism described here, can reduce 
convective precipitation, inducing and perpetuating higher probabilities of 
drought. The drought perpetuating influence is a recent development. In the 
past, the grazing animals would die off or be slaughtered in a drought. The 
pressure on vegetation would be reduced. Today support in fodder and water 
comes in a drought event, and vegetation continues to be under unrelenting 
impact, that is not alleviated by previous checks and balances. 
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Figure Captions 

Figure 1 The solar zenith angle 0 Q (the top graph), the absorption rate by 
the soil and the absorption rate B by the plants (in the middle), the 
enhancement ratio E_ of absorption by plants/soil surface to that by bare soil 

d 

and ratio P - B /(B + B ) (the lower graphs), all plotted vs time of day from 

P P 9 

sunrise (600h) to noon (1200h) at equinox, at the Sinai/Negev boundary. 
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TABLE 1. The soil absorptivity V the plant absorptivity b dp' the plant 
absorptivity as a fraction of the surface absorptivity, the enhancement 
factor E^ a by albedo reduction, the enhancement factor E dg by soil heat flux 
reduction (for g Q of 0.5 and 0.4) and the combined enhancement factor = 
E da* E dg ^ or ^o 0*4), all calculated for the direct solar beam 
only at different solar zenith angles. 


e Q ( degrees) 


30 

40 

45 

50 

60 

70 

b dg 


0.517 

0.490 

0.475 

0.457 

0.410 

0.335 

b dp 


0.194 

0.235 

0.259 

0.287 

0.360 

0.478 

b d 


0.711 

0.726 

0.734 

0.744 

0.770 

0.813 

P d 


0.273 

0.324 

0.353 

0.386 

0.468 

0.588 

b, +b, 

Ej 

da 0.58 


1.23 

1.25 

1.27 

1.28 

1.33 

1.40 

E *l + 9 ° Pd 

g o “ 0 -5 

1.27 

1.32 

1.35 

1.39 

1.47 

1.59 

d9 l-9 0 

V 0 * 4 

1.18 

1.22 

1.24 

1.26 

1.31 

1.39 

E da E dg 

g o =0.5 

1.56 

1.66 

1.71 

1.78 

1.95 

2.22 


g -0.4 

1.45 

1.52 

1.56 

1.61 

1.74 

1.95 
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TABLE 2. The absorption rate by the soil, the absorption rate Bp by the 

plants, the ratio P =■ Bp/(Bp+Bg) , the enhancement factor by albedo 

reduction E =(b +b )/0.58, the enhancement factor E by soil heat flux 
a P 9 9 

reduction (for g of 0.5 and 0.4) and the combined enhancement factor E=EiE„ 
3 o a g 

(for g Q of 0.5 and 0.4), all calculated for the global irradiance at different 
solar zenith angles. 


0 ( degrees) 



30 

40 

45 

50 

60 

70 

B 

g 



0.3745 

0.3081 

0.2716 

0.2333 

0.1544 

0.0786 

B p 



0.1454 

0.1505 

0.1498 

0.1466 

0.1314 

0.1004 

p 



0.280 

0.328 

0.356 

0.386 

0.460 

0.561 

E 

a 



1.23 

1.25 

1.27 

1.28 

1.32 

1.38 



/ g o = 0 -5 

1.28 

1.33 

1.36 

1.39 

1.46 

1.56 

9 0 P 

E -1+ 

9 1+ 9 0 

1 

! 

1 

V V 0 - 4 

1.19 

1.22 

1.24 

1.26 

1.31 

1.37 



r g o =0,5 

1.57 

1.66 

1.72 

1.78 

1.93 

2.16 

E=E 4 E g 

I 

) 

i 

I 

vg o =0.4 

1.46 

1.57 

1.61 

1.61 

1.73 

1.90 
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